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Abstract: We examine the possibility of a light (below 46 GeV) neutralino dark matter 
(DM) candidate within the 19-parameter phenomenological Minimal Supersymmetric Stan- 
dard Model (pMSSM) in the light of various recent experimental results, especially from 
the LHC, XENONIOO, and Planck. We also study the extent of electroweak fine-tuning for 
such a light neutralino scenario in view of the null results from the searches for supersym- 
metry so far. Using a Markov Chain Monte Carlo likelihood analysis of the full pMSSM 
parameter space, we find that a neutralino DM with mass > 10 GeV can in principle still 
satisfy all the existing constraints. Our light neutralino solutions can be broadly divided 
into two regions: (i) The solutions in the 10 - 30 GeV neutralino mass range are highly 
fine-tuned and require the existence of light selectrons (below 100 GeV) in order to satisfy 
the observed DM relic density. We note that these are not yet conclusively ruled out by the 
existing LEP/LHC results, and a dedicated analysis valid for a non-unified gaugino mass 
spectrum is required to exclude this possibility, (ii) The solutions with low fine-tuning are 
mainly in the 30 - 46 GeV neutralino mass range. However, a major portion of it is already 
ruled out by the latest XENONIOO upper limits on its spin-independent direct detection 
cross section, and the rest of the allowed points are within the XENONIT projected limit. 
Thus, we show that the allowed MSSM parameter space for a light neutralino DM below 
the LEP limit of 46 GeV, possible in supersymmetric models without gaugino mass unifi- 
cation, could be completely accessible in near future. This might be useful in view of the 
recent claims for positive hints of a DM signal in some direct detection experiments. 
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1 Introduction 

Low-scale Supersymmetry (SUSY) (see e.g., [1, 2]) is one of the most attractive candidates 
for New Physics beyond the Standard Model (SM). Apart from providing successful gauge 
coupling unification and a solution to the gauge hierarchy problem, it offers a natural 
candidate for Dark Matter (DM) in our Universe in the form of the lightest supersymmet- 
ric particle (LSP). In the -R-parity conserving Minimal Supersymmetric Standard Model 
(MSSM), the lightest neutralino (x?) is one of the most viable Weakly Interacting Massive 
Particle (WIMP) DM candidates (for a review, see e.g., [3, 4]). It can explain the observed 
DM relic density, while predicting experimentally accessible direct and indirect detection 
rates, over a wide range of supersymmetric model parameters, some of which are already 
getting constrained from the ongoing direct searches for the supersymmetric particles at 
the Large Hadron Collider (LHC), in combination with other low-energy data (for a review, 
see e.g., [5]). 

The recent hints of positive signals in three DM direct detection experiments, namely, 
DAMA [6, 7], CoGeNT [8, 9] and CRESST-II [10], have generated a lot of interest in light 
WIMP candidates in the 5-50 GeV mass range. This interpretation is however chal- 
lenged by the null results from various other direct detection experiments, most notably 
the latest XENONIOO results [11], which provide the most stringent upper limits on the 
spin-independent WIMP-nucleon scattering cross section for m^^ > 7 GeV ^. Nonetheless, 
due to the relatively poor sensitivity of the XENONIOO experiment in the very low WIMP 
mass regime, it is believed that an agreement between the positive and null sets of exper- 
imental results could be possible, if at all, only in this low mass region. Hence, it might 
be worthwhile examining the allowed MSSM parameter space to see if there exists a lower 
bound on the lightest neutralino mass irrespective of the direct detection results. 

^More recently, the TEXONO experiment [12] has achieved slightly better sensitivity than XENONIOO 
for m^<7 GeV. 
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The neutralino mass eigenstates in the MSSM result from mixing of the neutral 
bino (B), wino {W^) and higgsinos {H^,H^). This mixing is determined by the MSSM 
tan/3 parameter and the bino, wino and higgsino mass parameters Mi,M2 and respec- 
tively. In the SUSY models with gaugino mass unification at the Grand Unified Theory 
(GUT) scale, a relation between the bino and wino mass follows at the electroweak scale: 
Ml = I tan^ ewM2 ^ O.5M2 (see e.g., [13]) which, after mixing, translates into a chargino- 
neutralino mass relation. Therefore, a lower limit on the lightest neutralino (x?) mass of 
about 46 GeV can be derived for these models from the Large Electron Positron (LEP) 
chargino mass limit [14], whereas in the constrained MSSM (cMSSM) [15] with both gaug- 
ino and sfermion mass unification, this limit increases to well above 100 GeV from the 
strong constraints set by the recent LHC data [5]. 

On the other hand, in a generic MSSM scenario without the assumption of gaugino 
mass unification, there is no general lower limit on the lightest neutralino mass [16]. The 
LEP limit on the invisible decay width of the SM Z boson applies to light neutralinos 
with a mass below mz/2 = 45.6 GeV, but it depends on the ZxiXi coupling which could 
be small or even zero, depending on the higgsino component of the neutralino. In such 
a case, light neutralinos are mainly constrained by the DM relic density measurement as 
well as by the collider and flavor constraints on the SUSY parameter space. Assuming that 
the lightest neutralino is non-relativistic and provides the entire cold DM content of the 
Universe, while satisfying the LEP bounds on chargino and stau masses, Ref. [17] obtained 
a lower bound of nir-o > 18 GeV. This was relaxed to about 6 GeV without violating the 
LEP bounds and flavor sector constraints in SUSY models with a pseudoscalar Higgs boson 
(A) mass niA < 200 GeV and a large tan /3 [18, 19]. This was even further lowered to about 
3 GeV by allowing explicit CP violation in the MSSM Higgs sector [20]. 

Meanwhile, several new experimental results have been obtained at the LHC: (i) A new 
Higgs-like neutral scalar particle has been discovered [21, 22] with mass around 125 GeV 
which falls squarely within the MSSM-predicted range for the lightest CP-even neutral 
Higgs mass: rrih S [115,135] GeV (see e.g., [23]); (ii) The rare decay — )• jjL'^ ji^ was 
observed with a branching fraction of B{Bg — )• fi~^fi^) = (ji.2'^\'2) x 10^^ [24] which is 
in agreement with the SM expectation, B{B'^ A'"^M")sm = (3.23 ± 0.27) x 10"^ [25]; 
(iii) The lack of a SUSY signal at the ^/s = 7 and 8 TeV LHC has pushed the lower 
limits on the squark and gluino masses to about 1 TeV and beyond [26, 27]; (iv) Updated 
bounds have been obtained for the MSSM Higgs sector [28-30]. All these new results 
have profound implications for a light neutralino scenario within the MSSM, and some of 
these aspects have already been investigated in a number of recent studies [31-44]. The 
general conclusion is that light neutralino DM candidates with mass below about 15 GeV 
are severely constrained in generic MSSM scenarios (without gaugino mass unification). 

Another important issue to be addressed in the light of the recent LHC results is the 
apparent "little hierarchy problem", i.e., how does a multi-TeV SUSY particle spectrum 
conspire to give a weak-scale Z-boson mass and also a Higgs boson mass of 125 GeV? One 
way of analyzing this issue quantitatively is to evaluate the measure of electroweak fine- 
tuning (EWFT) by examining the minimization condition in the MSSM Higgs potential 
which determines the Z-boson mass [45, 46]. It is well-known that radiative corrections 
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must play a crucial role in determining the allowed SUSY parameter space necessary to 
generate a 125 GeV Higgs boson mass, much larger than its tree- level prediction of mh < 
mz- This in general can lead to a large fine-tuning. In addition to this, the requirement of 
a light neutralino DM will pose a challenge for any MSSM scenario, the severity of which 
is however strongly model-dependent. The naturalness of various SUSY models with a 
neutralino LSP has been analyzed in the literature (for an incomplete list, see [41, 47-60], 
and references therein). 

In this paper, we perform a dedicated study focusing on the naturalness of a light 
neutralino, and also examining how light the neutralino could be, after taking into account 
all the existing theoretical and experimental constraints. To perform such an analysis in 
the full 124-parameter MSSM is quite unrealistic, and hence, we need to make some well- 
motivated assumptions in order to reduce the number of parameters to a manageable level. 
Most of the earlier studies on SUSY focused on the cMSSM having only 5 parameters, 
assuming certain boundary conditions at the GUT-scale. However, in view of the latest 
null results from SUSY searches at the ^/s = 8 TeV LHC, the cMSSM seems too restrictive 
for low-scale SUSY phenomenology as the allowed cMSSM parameter space accessible to 
the ^/s = 14 TeV LHC is rapidly shrinking (for the latest global status, see e.g., Ref. [GI- 
GS]). Therefore, in this paper we choose not to make any assumptions at the high scale, 
and focus only on the low-scale MSSM parameter space from a phenomenological point of 
view. More precisely, we consider a CP-conserving MSSM (i.e., with no new CP phases) 
with Minimal Flavor Violation [G4] and with first two generations of sfermions degenerate. 
This is widely known as the phenomenological MSSM (pMSSM) [G5] (also known as 'SUSY 
without prejudice' [GG]) with 19 free parameters at the electroweak scale. We also study 
the level of fine-tuning for the light neutralino scenario in this context. 

In order to efficiently explore the 19-dimensional pMSSM parameter space, we perform 
a Markov Chain Monte Carlo (MCMC) likelihood analysis (for a review, see e.g., [G7]), 
with the priors optimized for a light neutralino scenario with mass below the conservative 
LEP lower limit of 46 GeV. We include in our analysis the latest experimental results for 
SUSY searches from the LHC [26, 27] which now supersede the Tevatron results [G8, 69], 
in addition to the existing LEP limits [14], wherever applicable. We also include the 
latest astrophysical/cosmological constraints for a WIMP DM from the 9-year WMAP 
data [70] as well as the very recently released Planck data [71]. We further examine the 
allowed parameter space in the light of various recent results for DM direct detection, most 
notably the XENONIOO limits [11], as well as the indirect detection results from Fermi-LAT 
data [72]. 

We find that a light neutralino DM with mass as low as 10 GeV is still allowed in 
the pMSSM, while satisfying all the existing experimental constraints provided we only 
take the model-independent analysis results from LEP. (Including the LEP limits strictly 
applicable to gaugino-mass unification models allows only the solutions with m~o > 30 
GeV, in agreement with previous results [33, 34].) However, such neutralinos which are 
required to be mostly bino-like are severely fine-tuned and require the existence of light 
sleptons below 100 GeV in order to provide an efficient annihilation channel to reduce 
the bino relic density to be consistent with the observed limit. A dedicated analysis of 
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the LEP data in the context of a pMSSM scenario could ehminate this region completely. 
On the other hand, low fine-tuning regions can be obtained around m^o = 45 GeV where 
the resonant annihilation via the s-channel Z-exchange is possible for the neutralino with 
a non-negligible higgsino component. However, such regions are mostly excluded by the 
recent XENONIOO limits on the spin-independent DM-nucleon scattering cross section, and 
the remaining such points are well within the reach of the XENONIT [73] and LUX [74] 
projected limits. 

Our paper is organized as follows: In Section 2, we briefly discuss the electroweak fine- 
tuning measure. In Section 3, we list all the pMSSM parameters and their scan ranges, and 
also summarize all the relevant experimental constraints used in our numerical analysis. In 
Section 4, we present our scan results and discuss their implications for a light neutralino 
DM. Finally, our conclusions are given in Section 5. 



2 Electroweak Fine- Tuning 



To quantify the amount of fine-tuning in the electroweak symmetry breaking (EWSB) 
sector of the MSSM, it is sufficient for us to analyze the tree-level MSSM scalar potential. 
Since it allows only charge-conserving vacua, we only have to minimize the scalar potential 
for the neutral scalar fields [1]: 
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where fi is the SUSY-preserving bilinear Higgs superpotential parameter, ^ and B are 
the soft scalar masses and the bilinear coupling in the SUSY-breaking sector respectively, 
and g,g' the SU{2)l and U{1)y gauge couplings, respectively. After minimization, we 
obtain the well-known relation for the Z-boson mass: 
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where tan/3 = Vu/vd is the ratio of the vacuum expectation values (vevs) of the two Higgs 
doublet fields Hu and respectively. It is clear from Eq. (2.2) that a cancellation of the 
terms on the right hand side (RHS) is required in order to obtain the measured value of 
mz = 91.2 GeV [5] especially if the mass parameters on the RHS are orders of magnitude 
larger than the weak scale which indeed seems to be the case, given the current experimental 
status of the direct SUSY searches [26, 27]. Thus naively speaking, the weak scale value of 
can be used as a measure of fine-tuning in the MSSM. 

A more sophisticated way to quantify the degree of EWFT is by using log-derivatives [46] : 
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where piS are the parameters that determine the observable Z-mass at tree- level. From 
Eq. (2.2), we have pi = {fj? ,b,mH^,mHa} (with b = Bfi), and the total measure of the 
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EWFT is defined as 
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(2.4) 



with the individual Apj's given by [75] 
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Here m\ = i?/x(tan/3+cot /3) is the MSSM pseudo-scalar Higgs mass, and we have assumed 
tan /3 > 0. Values of Atot ^ 1 indicate significant fine-tuning. Note that in the limit of 
large tan/3, the quantities Amjj^ and Am?^^ in Eq. (2.5) are small, and 
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Thus in the limit of large tan /3, we recover the naive result that fine-tuning increases with 
increasing 

We should note here that including loop corrections to Eq. (2.2), one finds the largest 
contribution to be coming from the (s)top loop which feeds into the soft mass of the 
Higgs [13]: 
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where yt is the top-quark Yukawa coupling, At is the third-generation j4-term in the SUSY- 
breaking sector, and A is some high scale where the stop masses m-jr^ ^ are generated from 
the SUSY-breaking mechanism. Even for a low-scale SUSY-breaking scenario (such as 
gauge mediation), this requires a fine-tuning of at least a few percent in order to get the 
observed Z-mass and it becomes worse for heavy stop masses [76]. A heavy physical stop 
mass > 500 GeV is anyway required to provide large enough radiative corrections to 
the light CP-even Higgs mass to raise it to the vicinity of 125 GeV from its tree-level value 
<mz [13]: 
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where Xt = At — ^ cot /3 is the stop mixing parameter. A fine-tuning measure A^^ for the 
Higgs mass can be defined analogous to Eq. (2.3), and for a particular choice of some of the 
SUSY parameters, it was found that A^^ > 75 (100) in order to achieve a Higgs mass of 
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124 (126) GeV and the corresponding stop mass is always heavier than 300 (500) GeV [77]. 
The fine-tuning due to Eqs. (2.7) and (2.8) could in principle be added in quadratures to 
our "total" fine-tuning parameter defined by Eq. (2.4), but since their specific values are 
scale-dependent ^ and adds to some arbitrariness in its definition, we do not include them 
in our analysis. 

3 Parameters and Constraints 

We consider the pMSSM with 19 free parameters at the SUSY scale, as shown in Table 1. 
We further assume that the lightest neutralino is the LSP, and our goal is to examine 
how light and natural the neutralino could be while satisfying all the existing theoretical 
and experimental constraints on the MSSM parameter space. In particular, we focus 
on the lightest neutralino masses below the LEP limit of 46 GeV [14] which is strictly 
valid assuming gaugino mass unification at the GUT-scale. We perform a numerical scan 
over the 19-dimensional pMSSM parameter space using a MCMC-based likelihood analysis 
optimized for such light neutralino scenarios with the prior ranges given in Table 1. These 



Parameter 


Description 


Prior Range 


tan/? 


Ratio of the scalar doublet vevs 


[1, 60] 




Higgs-Higgsino mass parameter 


[-3, 3] TeV 


Ma 


Pseudo-scalar Higgs mass 


[0.3, 3] TeV 


Ml 


Bino mass 


[-0.5, 0.5] TeV 


M2 


Wino mass 


[-1, 1] TeV 


M3 


Gluino mass 


[0.8, 3] TeV 




First/second generation Ql squark 


[0, 3] TeV 




First /second generation Ur squark 


[0, 3] TeV 




First/second generation Dji squark 


[0, 3] TeV 




First /second generation Ll slepton 


[0, 3] TeV 


men 


First/second generation Eji slepton 


[0, 3] TeV 




Third generation Ql squark 


[0, 3] TeV 


m-r 


Third generation Uji squark 


[0, 3] TeV 


mr 

OR 


Third generation D/j squark 


[0, 3] TeV 




Third generation Ll slepton 


[0, 3] TeV 




Third generation Eji slepton 


[0, 3] TeV 


At 


Trilinear coupling for top quark 


[-10, 10] TeV 


Ab 


Trilinear coupling for bottom quark 


[-10, 10] TeV 


Ar 


Trilinear coupling for r-lepton 


[-10, 10] TeV 



Table 1. The pMSSM parameters and their range of values used in our numerical analysis. 

particular ranges are chosen in order to economize the scan time and to focus only on the 
^The SUSY-scale by convention is usually taken to be the geometric mean of the two stop masses. 
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SUSY parameter space not yet disfavored by the combined direct search results from the 
LEP [14], Tevatron [68, 69] and LHC [26, 27], as discussed later in this section. For the 
SM parameters as{mz), aem("iz), mw, fnt and m^, we use the standard values as given 
in Ref. [5]. 

Apart from the direct collider constraints on the sparticle masses, there exist various 
theoretical and experimental constraints which must be imposed on the pMSSM parameter 
space in our analysis. As a standard theoretical requirement, our sparticle spectrum for 
each allowed point in the parameter space must be tachyon-free and should not lead to color- 
and charge-breaking minima in the scalar potential [78]. We also require that the scalar 
potential is bounded from below and is consistent with electroweak symmetry breaking. 
From the radiative electroweak symmetry breaking arguments, we can restrict the tan /3 
parameter to be roughly between 1 - 60, as given in Table 1. 

The various experimental constraints from direct collider searches, Higgs and flavor 
sectors, and astrophysical/cosmological data used in our analysis are summarized below. 

Invisible Z-decay Width 

The precise measurement of the Z-boson decay width at LEP: F^* = 2495.2 it 2.3 MeV [79] 
puts severe constraints on light neutralinos and charginos with mass < mz/'2- Prom the 
LEP measurements of the invisible decay width of the Z-boson: F^"^ = (499.0 it 1.5) 
MeV [79] , the parameter space for the lightest neutralino in our case is restricted to mostly 
gaugino-like scenarios (|/i| ^ Mi^2) since the neutralino coupling to Z is only via its higgsino 
component. The allowed fraction of the higgsino component for a given neutralino mass 
can be calculated using the following expression for the partial decay width of the Z-boson 
to neutralinos [80]: 

F(Z ^ xgxg) 
F(Z ^ vv) 

where T[Z — )• vv) = (501.62 it 0.10) MeV is the SM contribution to its invisible decay 
width (for 3 neutrino species). The components O^-^ are defined as 

Ofj = + \n,,N*^ , Of^ = - {Of-Y (3.2) 

with Nij measuring the gaugino-Higgsino mixing: 

4 

X^ = Y,N,kipl, where V^o = (5, #0,^0). (3.3) 

fe=i 

Using the LEP measurement of the invisible Z-decay width, the following constraint 
can be derived: 

F(Z ^ X°Xi) < 3 MeV if (m^, + m^J < mz (3.4) 

Note that this constraint should apply to all light neutralinos satisfying this condition, 
not just the LSP (i = j = 1 in Eq. 3.4). However, it is unlikely that decays such as 

- 7- 
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Z — >■ X2X1 will be kinematically allowed, and in such cases, X2 "^i^ mostly decay to visible 
final states. Similarly, the decays Z — t- vivj are not kinematically allowed for the parameter 
space examined here, and hence, they do not contribute to the purely invisible width of 
the Z-boson. 

Exclusion Limits from Collider Searches 

The experimental lower limits on the sparticle masses are usually quoted assuming gaugino 
and/or sfermion mass universality at the GUT scale. In a generic MSSM setup, most of 
these constraints can be relaxed, or even circumvented, for example in case of small mass 
splitting with the LSP or in case of small couplings to the SM vector bosons. Since we 
are interested in light neutralino LSPs here, we must carefully interpret the direct search 
limits in order to be able to include all the allowed pMSSM parameter space. 

Neutralino: As already mentioned in Section 1, there are no rigorous lower limits on 
the neutralino masses in the MSSM from direct collider searches. The LEP limits [14, 81] 

m~o > 46 GeV, m~o > 62.4 GeV, m~o > 99.9 GeV, m~o > 116.0 GeV (3.5) 

were derived assuming gaugino mass unification at the GUT scale, and hence, relating the 
neutralino mass to the chargino mass. Moreover, for a (mostly) bino-like neutralino, which 
is required to be the case for m^o < in order to avoid the Z-width constraint, its 

production via s-channel exchange of Zj^* is (negligible) absent. The t-channel production 
cross-section via selectron exchange is also expected to be small for selectron masses above 
the LEP limit (see below). Thus, we can easily satisfy the LEP upper limits on the 
neutralino pair-production cross sections (T(e^e~ — )• X^Xj) [81, 82] for a mostly bino-like 
neutralino LSP. Similarly, the Tevatron [83] and LHC [84] SUSY searches for final states 
involving Z-bosons cannot constrain a bino-like neutralino. 

Chargino: Charginos can be pair-produced at LEP via s-channel exchange of Z/7* 
or f-channel exchange of electron-sneutrino, with destructive interference. It dominantly 
decays to G)^ if kinematically allowed. If not, the three-body decay to ff'Xi via virtual 
T^-boson or sfermions becomes important in which case, the final state fermions (/, /') 
are dominantly leptonic (hadronic) if the sleptons are light (heavy). From the combined 
searches in fully-hadronic, semi-leptonic and fully-leptonic decay modes, LEP has derived a 
general lower limit of 103.5 GeV [14] which is valid for pMSSM as well, except in corners of 
phase space where (i) the detection efficiencies are reduced, e.g., when the mass differences 
A7n+ = m~± — nir-o or Am^ = m~± — rriy are very small (below a few GeV); or (ii) the 
chargino production cross section is suppressed, e.g., when the electron sneutrino mass 
is small, thus leading to a destructive interference between s- and t-channel Feynman 
diagrams. Dedicated searches for such scenarios have also been performed. For instance, 
for small Am+ < 3 GeV but with large sneutrino mass, the limit becomes m~± > 91.9 

Al 

GeV for degeneracy in the gaugino region (|Mi| ~ IM2I ^ while m~± > 92.4 GeV for 
degeneracy in the higgsino region (|/i| <C IM2I) [14]. Without assuming gaugino mass 
unification, a lower limit of m~± > 70 GeV was set for any Am^ and m-p > 300 GeV. For 

Al 

smaller sneutrino masses, the sensitivity decreases due to the reduced pair production cross 
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section and also due to reduced selection efficiency. In such situations where none of the 
above mass limits can apply, the generic lower limit of approximately 45 GeV, derived from 
the analysis of the Z-width, is still valid since this is independent of the field composition 
and of the decay modes of the charginos. Note that unlike neutralinos which couple to 
the Z-boson only via their higgsino component, the charginos couple to Z via both their 
gaugino as well as higgsino components; so it is not possible to avoid the Z-width constraint 
for a light chargino. 

Sneutrino: Light sneutrinos can only decay invisibly to z^Xi unless the decays to 
charginos and heavier neutralinos are not kinematically suppressed. The invisible width 
of the Z-boson puts a lower limit on the left-sneutrino mass of 43.7 GeV, which improves 
slightly to 44.7 GeV if all three sneutrinos are mass-degenerate. Note that the lightest left- 
sneutrino by itself cannot be a cold DM candidate [85, 86], and we must introduce a mixing 
with a SM singlet sneutrino to make it a viable DM candidate (see, e.g., [87-93]). Since 
we are dealing here only with the MSSM field content and do not have a right-sneutrino 
component, we discard those points for which the sneutrino is the LSP. 

Slepton: Studies of the Z-boson width and decays put a lower bound on the slepton 
masses "^J^^^j > ^0 (41) GeV, independently of the decay modes for individual sleptons 

{t = e, /I, r). This limit improves to 43 GeV if all the three slepton flavors are mass- 
degenerate. Tighter limits can be obtained assuming that sleptons are pair-produced at 
LEP, and each slepton dominantly decays to ^x?) ^^^^s leading to two back-to-back leptons 
and missing transverse momentum. These limits are valid for a mass splitting Am^ = 
mj — rn:^o > 15 GeV so that the final state leptons are not too soft. Moreover, the LEP 
results are interpreted assuming that only £r£r production contributes, and hence, the 
limits are usually quoted for Ir, since it is typically lighter than in most SUSY models, 
and has a weaker coupling to the Z-boson so that the limits are more conservative. This 
is a good approximation for selectrons and smuons, but not for staus which can have 
significant mixing between the fiavor eigenstates tx and tr. The most conservative limit 
on the mass of the lightest stau is obtained with a mixing angle ~ 52° which minimizes 
the production cross section. 

The slepton mass limits of 0(100) GeV quoted in Ref. [14] were derived under the 
assumption of gaugino mass unification at the GUT scale which was used to fix the masses 
and composition of neutralinos. It was also assumed that the slepton branching ratio 
to ^Xi is nearly 100% which is a good approximation if the second lightest neutralino is 
heavy enough to suppress the cascade decay into £x2 followed by X2 ~^ ffXi or Xil- 
For smuons and staus, the LEP limits are independent of the MSSM parameters [94], 
and hence, directly applicable to our case. However, the selectron mass limit will be 
different if we do not assume gaugino mass unification due to a different production cross 
section involving the t-channel neutralino exchange, in addition to the usual s-channel 
7*/Z exchange. Since there is no dedicated analysis of the LEP data addressing this issue 
available in the literature and this is beyond the scope of our present work, we only use 
the generic lower bound for e derived from the Z- width in our numerical analysis, but will 
also comment on the implications of the tighter selectron mass bound from LEP on our 
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light neutralino scenario. We will also include the latest 95% C.L. LHC exclusion limits for 
slepton pair production interpreted in the slepton-neutralino mass plane of pMSSM [95] 
which are applicable to both left- and right-handed selectrons and smuons. 

Gluino and Squarks: The current LEP limits on the squark masses are similar to 
the slepton mass limits of 0(100) GeV. However, since squarks are colored objects, their 
production cross sections are much higher at hadron colliders. The highest sensitivity on 
squark and gluino production now comes from the LHC experiments. The generic lower 
limit on the first / second generation squark masses is 600 - 750 GeV and on the gluino mass 
is 700 - 900 GeV, as set for simplified SUSY models by the ATLAS analysis of the 8 TeV 
LHC data [26]. The corresponding CMS limits are very similar [27]. However, from the 
latest global fit of pMSSM after the LHC results [59] (see also Ref. [96, 97] for an earlier 
global fit of pMSSM), the corresponding lower limit on squark mass is rrig > 500 GeV, and 
hence, we use this value to constrain our pMSSM parameter space. 

The LHC and Tevatron limits on the third generation squarks are usually weaker since 
the amount of bottom and top quark content in the proton is negligible, and hence, the 
direct production of bottom and top squark is suppressed with respect to the first/second 
generation squarks. The current exclusion limit for top squarks is between 300 - 600 GeV 
from the LHC data, depending on the decay channel [26, 27]. Similar limits have been 
derived for the bottom squark as well. Following the latest global fit of pMSSM [59], we 
set the lower limits for light top squark and sbottom masses at 400 GeV and 300 GeV, 
respectively. Note that these limits are applicable as long as the light stop/sbottom is 
not highly mass degenerate with the lightest neutralino (with Am < 10 GeV or so) which 
always turns out to be the case for our light neutralino solutions satisfying all the other 
constraints. Thus we do not have any light sbottom solutions as considered in Ref. [38]. 

The lower limits on the sparticle masses derived from the above discussion are summa- 
rized in Table 2. We emphasize here that for a light neutralino DM with m~o < mz/2 as 
considered in our case, the dominant annihilation channels will be the t-channel processes 
mediated by light sfermions. The neutralino/chargino co-annihilation is also important if 
Am+ is small. Hence, the lower limits on the chargino and sfermion masses as given in 
Table 2 are crucial ingredients in our numerical analysis. 

Bs — >■ 

The rare decay Bg — s- fJ-^ fJ^^ is known to be a sensitive channel for new physics since its 
SM predicted rate is small due to helicity suppression. Recently, the first evidence for 
this decay was observed by the LHCb collaboration with the measured branching ratio: 
B{Bs fJ-^fJ-') = i^-'^-li) X 10"^ [24] which is in excellent agreement with the SM 
prediction of (3.23 ± 0.27) x 10"^ [25], thus raising some concerns for the "health" of 
SUSY. However, it must be noted that the upper limit derived from the latest LHCb result 
is actually slightly weaker than the earlier upper limit of < 4.5 x 10~^ [98]. The effect 
of the new results is mostly felt in the large tan /? > 50 regions of the MSSM parameter 
space [99] which are also strongly constrained from direct SUSY searches and the MSSM 
Higgs searches. 
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Particle 


Mass limit (GeV) 


Validity Condition 


xt 


103.5 


m~+ — rrizo > 3 GeV, mj > m^± 




70 


m~ > 300 GeV l/vl > \Mo\ 




45 


generic LEP bound 


H'K 


88 


rr)~ _ jTn^f, > 1 5 CipV fTT — >■ = 1 






rri- — m-n ^ 1 ^ HpV RR ^t^ —V tv^^ — 1 

iit'j--^ //6~0 ^ XO Vjc V , 1 J iV / \ t ' — ^ 


e-R 


yo 








= — zUU GeV, tanp = 1 


^R (L) 


40 (41) 


crp'nprir' TiT^/P* noiinrl 


V 


43.7 


generic LEP bound 


9 


800 




q 


500 




h 


400 




bi 


300 





Table 2. The lower limits on the sparticle masses used in our numerical analysis. The chargino and 
slepton mass limits are derived from the LEP data [14] while the squark and gluino mass limits are 
derived from the LHC data [26, 27] which now supersede the LEP as well as the Tevatron [68, 69] 
limits. 

6 — >■ S7 

We have also included the constraint from the radiative i?-meson decay branching ratio, 
B{b 57) = (3.55 ± 0.24 ± 0.09) x 10""^ [100] which is somewhat higher than the SM 
prediction of (3.15 it 0.23) x 10~^ [101]. Thus large SUSY corrections are preferred which 
mainly occur for light chargino and top squarks and for large tan /3 [102] . 

Another important constraint comes from the muon anomalous magnetic moment = (g— 
2)^/2 which gives a more than 3a discrepancy with the SM prediction: = (26.1 ±8.0) x 
10-1° [103]. The SUSY contribution to da^ can explain this discrepancy with relatively 
light smuons and/or large tan/3 [104-106], and < region of the SUSY parameter space 
is strongly disfavored (unless the electroweak gaugino masses Mi, M2 < 0). 

Higgs Sector 

For the light CP-even Higgs mass, we have chosen the value of ruh = (125 it 2) GeV, 
following the latest best-fit mass measurements of the Higgs-like particle discovered at the 
LHC: (125.2 ±0.3 ±0.6) GeV (ATLAS) [107] and (125.8 ± 0.4 ± 0.4) GeV (CMS) [108]. 
For the other MSSM Higgs bosons, we ensure that all our allowed points satisfy the latest 
LHC constraints on the — tan/3 plane [109, 110] and on the MSSM charged Higgs 
mass [111, 112] which are related to each other at tree-level by rnjj^ = m^+my^. Note that 
the non-decoupling region with light itla ~ 95 — 130 GeV, almost mass-degenerate with the 
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light CP-even Higgs, and with the heavy CP-even Higgs SM-Uke, is disfavored [113, 114] by 
the latest LHC Higgs data and flavor constraints, especially the Bg — )• fJ''^f^~ and b — >• 57. 
Hence, we work in the decoupling region with niA > 300 GeV with the light CP-even Higgs 
SM-like with mass in the vicinity of 125 GeV, and with the heavy CP-even Higgs nearly 
mass-degenerate with the CP-odd Higgs. Also note that for a mostly bino-like neutralino, 
the hxiXi coupling is small enough to satisfy the current global limit on the Higgs invisible 
decay branching ratio: BRinv < 0.28 at 95% CL [115]. 

Dark Matter Constraints 

The latest results from Planck give the current relic density of the cold dark matter content 
m our universe to be n^h^ = 0.1199 ± 0.0027 at 68% CL [71]. The corresponding value 
from the 9-year WMAP data is Q^h"^ = 0.1148 ± 0.0019 [70]. For the relic density of 
the neutralino DM in our case, we only require it to satisfy the WMAP 2a upper bound 
combined with 10% theoretical uncertainty: ^^oh'^ < 0.138 which also encompasses the 
latest observed value from Planck. The cases where the neutralino relic density is below 
the corresponding WMAP lower bound, il:^oh? < 0.091, could correspond to some multi- 
component DM scenario (see e.g., [116-125]). 

As far as the DM direct detection constraints are concerned, since there is no unani- 
mous upper bound on the direct detection cross section, we do not put this constraint a 
priori on the model parameter space. However, as we will see later, most of the allowed 
parameter space satisfying the other constraints also satisfy the most stringent upper limit 
on the spin-independent DM-nucleon scattering cross section set by the XENONIOO ex- 
periment [11]. Note that for > 10 GeV, the spin-independent direct detection con- 
straints are more stringent than the collider constraints from monojet [126, 127] and iso- 
lated monophoton CMS-monoph, ATLAS-monoph searches at the LHC as well as from the 
LEP mono-photon data [130]. 

Complementary to the direct detection constraints, there exist indirect detection con- 
straints which are mostly sensitive to light WIMPs annihilating to SM fermions which 
eventually lead to gamma-ray signals. A lower limit of m~o > 10 GeV was derived from 
the CMB constraints [131-135]'^, assuming a velocity-independent annihilation cross section 
of {a av) — 3 X 10 ^'^cm"^ • s ^. Under the same assumption, the Fermi-LAT data put lower 
bounds of m~o > 27 GeV for annihilation to bb channel and nir-o > 37 GeV for t'^t~ chan- 
nel [137]. However, these bounds can be relaxed if we include the velocity-dependent contri- 
butions, as shown for pMSSM in Ref. [138]. In our numerical analysis, we include the latest 
Fermi-LAT 95% CL upper limit on the integrated 7-ray flux: (p^ < Ax lO^^'^cm"^ -s"^ [72]. 

The various experimental constraints discussed above and used in our numerical anal- 
ysis, in addition to the sparticle mass limits listed in Table 2, are summarized in Table 3. 



^Similar limits have also been derived from BBN constraints [136]. 



- 12 - 



Parameter 


Constraint 


■pinvisiblc 
^ Z 

Bib s-f) 
5a ^ 


(125 ± 2) GeV 
< 3 MeV 
< 0.138 
< 4 X 10"^° cm-2 • s"^ 
(3.2^1:^) X 10-9 
(3.55 ±0.26) X 10"^ 
(26.1 ±8.0) X 10-1° 



Table 3. The relevant experimental constraints used in our analysis, in addition to those listed in 
Table 2. 

4 Results 

In order to scan the 19-dimensional pMSSM parameter space more efficiently while sat- 
isfying all the constraints listed in Tables 2 and 3, we have performed a MCMC analysis 
using a Gaussian distribution of likelihood function: /(x, xq, cj) = exp[— (x — xo)^/2cr^] for 
all the observables, with a preferred value xq ± cj. We have used CaclHEP2.3 [139] and 
micrOMEGAs2.4 [140-142] to compute ah the observables, together with SoftSUSY [143] for 
calculating the particle spectrum. 

First we discuss our MCMC scan results for the relic density of a light neutralino 
DM candidate as shown in Figure 1. We require all the allowed points (shown as circles) 
to satisfy the experimental constraints given in Table 3, along with the LEP limits on 
sparticle masses given in Table 2. The latest LHC results put much tighter bounds on 
the strongly interacting squarks and gluinos and further eliminate some of these otherwise 
allowed parameter space, as shown by the starred points in Figure 1. The WMAP-9 2(T 
band is shown in grey, whereas the latest Planck result is shown as dark shaded region. We 
find that light neutralinos with mass as low as 10 GeV are still allowed, though severely 
fine-tuned with the electroweak fine-tuning measure defined by Eq. (2.4): Atot ^ 1- This 
can be understood as follows by analyzing the gaugino and higgsino components of the 
lightest neutralino as well as its dominant annihilation channels. 

The bino, wino and higgsino fractions of the lightest neutralino for all the allowed 
points in our pMSSM parameter space are shown in Figure 2. We reproduce the well- 
known result that the lightest neutralino is mostly bino-like for masses below m^/2, mainly 
due to the invisible Z-decay width constraint in Eq. (3.4). However, a purely bino DM 
tends to overdose the universe since its annihilation cross section is p-wave suppressed 
and is typically much smaller than the required thermal WIMP annihilation rate of 3 x 
IQ-^^ cm^ • One alternative is to consider a "well-tempered" neutralino [144] which 
corresponds to the boundary between a pure bino and a pure higgsino or wino. Another 
possibility to reduce the bino relic density is by annihilation via the t-channel slepton 
exchange (the so-called "bulk region") which is efficient for light sleptons, or by using co- 
annihilation with a light stau, stop, sbottom, chargino or second-lightest neutralino. We 
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Figure 1. The relic density of a light neutralino DM in pMSSM satisfying all the experimental 
constraints discussed in Section 3. The color-coding denotes the fine-tuning measure defined by 
Eq. (2.4). The points denoted by circles satisfy all the experimental constraints, except that the 
squark masses are only required to satisfy the LEP lower limits. For the starred points (a subset 
of the circled points), the corresponding squark masses satisfy the latest LHC constraints. The top 
(bottom) grey horizontal line shows the 2cr upper (lower) limit of the cold dark matter relic density 
from WMAP-9 data, whereas the black (shaded) region shows the la allowed range from the recent 
Planck data. 

find that most of the points with mc;o close to 45 GeV can have a non-neghgible higgsino 

'''-1 

component, and hence, can easily satisfy the WMAP upper limit on the rehc density since 
there is a resonant annihilation via Z-mediated s-channel. These points are also less fine- 
tuned. On the other hand, the light neutralino DM points in the 10 - 30 GeV range as 
shown in Figure 1 have to be mostly bino-like and lie either in the bulk or co-annihilation 
region, thus leading to significant fine-tuning. 

This is further clarified in Figure 3 where we show the various next-to-lightest super- 
symmetric particles (NLSPs) and their masses as a function of the lightest neutralino mass. 
We see that all the allowed points with rrir^o < 30 GeV have a charged slepton NLSP with 
mass below 100 GeV. Especially the points with a light stau are severely fine-tuned since 
they usually require a mass suppression by the off-diagonal elements in the slepton mass 
matrix, or a large //-term. We also show in Figure 3 the LEP exclusion regions in the 
charged slepton-neutralino mass plane, derived under the assumption of gaugino mass uni- 
fication [14]. The limits for light smuons and staus are still applicable to the pMSSM case 
as long as Ame > 15 GeV, but not directly to light selections if we assume non- universal 
gaugino masses, and hence, can still allow the low neutralino mass regime. The latest 
95% C.L. ATLAS exclusion limits [95] are also shown in Figure 3 which were derived from 
searches for direct slepton (selection and smuon) pair production and interpreted in the 
pMSSM. A similar dedicated analysis of the LEP data is required in order to completely 
rule out the light selectrons, and hence, the lightest neutralino DM mass below 30 GeV for 
the pMSSM scenario. 

The neutralino DM-nucleon spin-independent scattering cross-sections for the allowed 
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Figure 3. The various NLSP masses as a function of the LSP mass for the allowed points (circles) 
shown in Figure 1. The LEP exclusion regions strictly applicable for Jlh (red shaded) and ti (green 
shaded) and the LHC exclusion region for £ (l — e, /i) (blue shaded) are also shown. 



points is shown in Figure 4. It is interesting to note that most of the allowed region in 
Figure 1 with low fine-tuning lead to a higher scattering cross section via Z-boson exchange 
and are already ruled out by the latest XENONIOO data [11]. The projected limit of 
XENONIT experiment [73] (and also LUX [74]) will be able to completely eliminate the 
low fine-tuning region for a light neutralino DM in pMSSM. In the light of the recent claims 
for positive hints of a light DM from some experiments, it is worth mentioning here that 
a few of our solutions with a^^ ~ 10^'' pb are in the vicinity of the 2a preferred range of 
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Figure 4. The spin-independent direct detection cross-section values for the aUowed points in 
our pMSSM scan. The cofor-coding and labeUng of the points are the same as in Figure 1. The 
current upper hmit from the XENONIOO experiment (upper soUd hne) and the projected hmit from 
the XENONIT experiment (fower sohd Hne) are shown as sohd hnes. The 2(7-preferred range of 
CRESST-II is shown as the shaded region. 

the CRESST-II results [10], but are ruled out by the XENONIOO data. 

Figure 5 shows the integrated photon flux from the neutralino DM annihilation. Note 
that since the LSP in our case is mostly bino-like with heavy squarks and higgsinos and 
the charginos are not mass-degenerate with the LSP, the annihilation to photons is loop- 
suppressed, and hence, the photon line emission will be small. It is clear from Figure 5 
that the current sensitivity of the Fermi-LAT [137] still leaves most of our allowed pa- 
rameter space untouched. The future data from ongoing Fermi-LAT and next generation 
gamma-ray searches might be able to probe our allowed parameter space with a photon 
line signal [145]. 

For completeness, we also show in Figure 6 some other relevant pMSSM parameters 
with respect to the lightest neutralino mass. The niA — tan f3 parameter space is mostly 
consistent with the latest MSSM Higgs sector limits from the LHC [109, 110]. As for the 
bino mass parameter Mi, it is clear that a relatively small value of \Mi\ < 100 GeV is 
preferred to obtain a light bino-like neutralino LSP. Finally, as is well-known, a relatively 
large value of l^tl is required in order to enhance the radiative corrections for the light 
CP-even Higgs mass to be consistent with the LHC-preferred value of 125 it 2 GeV. 

Finally, we wish to point out that the allowed sparticle spectra discussed here with 
mostly heavy squarks, and with light sleptons, chargino and bino-like neutralino LSP, are 
also crucial for explaining the muon {g—2) anomaly [146], while simultaneously satisfying all 
the other experimental constraints. Light staus with sizable left-right mixing can also lead 
to an enhanced /i — )• 77 decay rate [147-149], which might be able to explain the persistent 
excess in the Higgs signal strength in this channel: 1.65 ± 0.24(stat)lQ^g(syst) [150]. 



-16- 



10-' 

10-10 k 

10-12 



>- 

-e- 



10 



16 



10" 



Fermi-LAT 



« • • 



18 I- LEP • 
ILHC * 



10 15 20 25 30 35 40 45 
mjjo [GeV] 



1000 



100 o 
< 



10 



Figure 5. The integrated photon flux from annihilation of the neutraUno DM as a function of its 
mass. The color-coding and labeling of the points are the same as in Figure 1. The solid horizontal 
line shows the current upper limit from the Fermi-LAT data. 
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Figure 6. The allowed parameter space for some of the relevant pMSSM parameters consistent 
with a light neutralino DM. The color-coding and labeling of the points are the same as in Figure 1. 



5 Conclusion 

We have studied the naturalness of a light neutralino dark matter candidate in the MSSM in 
the light of the latest results from the collider (LEP, Tevatron and LHC), flavor (LHCb) and 
dark matter (XENONIOO, WMAP-9, Planck) sectors. In particular, keeping in mind the 
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recent positive hints for a light DM below 45 GeV in some direct detection experiments and 
the null results from SUSY searches at the LHC, we perform a dedicated study focusing 
on the possibility of a light DM candidate in the form of the lightest neutralino in the 
pMSSM, also analyzing the naturalness of such a scenario. We include all the new results 
from the ongoing LHC experiments in our analysis, such as the null results from the SUSY 
searches, the observation of a Higgs-like particle, the observation of the rare meson decay 
Bs — )• , and the updated constraints on the MSSM Higgs sector. We also take into 

account the 9-year WMAP results as well as the recently released Planck results for the 
DM relic density and the Fermi-LAT data for the integrated photon flux. In order to 
efficiently scan over the 19-dimensional pMSSM parameter space, we perform a MCMC 
likelihood analysis optimized for a light neutralino with mass below the LEp bound of 46 
GeV (applicable to gaugino mass unification models). 

We find that a light neutralino DM with mass as low as 10 GeV is still allowed in 
the pMSSM, while satisfying all the existing experimental constraints. However, such 
neutralinos which are required to be mostly bino-like are severely fine-tuned and require 
the existence of light sleptons with mass below 100 GeV in order to provide an efficient 
annihilation channel to reduce the DM relic density below the observed upper limit. Such 
light smuons and staus are excluded from LEP searches while light selectrons are excluded 
only if we assume gaugino mass unification. A dedicated analysis of the LEP data in the 
context of a pMSSM scenario could completely eliminate the possibility of a light neutralino 
DM in the mass range of 10 - 30 GeV. We also find that for the allowed parameter space, 
other possible solutions in this mass range as discussed in the literature (e.g., light sbottom 
NLSP) are now excluded mainly due to the latest LHC results on the strongly interacting 
sfermion sector. 

On the other hand, low fine-tuning regions can be obtained around mc;o = 45 GeV 
where the resonant annihilation via the s-channel .Z-exchange is possible for a bino-higgsino 
mixture of neutralino LSP. However, such regions also predict a higher spin-independent 
DM-nucleon scattering cross section, and are mostly excluded by the recent XENONIOO 
limits. The remaining such points are within the reach of the XENONIT projected limits. 

In conclusion, while a light neutralino DM with mass below 46 GeV has been conclu- 
sively ruled out in MSSM with gaugino mass unification by LEP searches, such a possibility 
in a general version of MSSM is of enormous interest in the light of the recent claims for 
a positive signal in some DM direct detection experiments. Taking into account the latest 
experimental results from collider, flavor, dark matter and astrophysical/cosmological sec- 
tors, we show that such a light neutralino DM scenario is also getting highly constrained in 
MSSM without gaugino mass unification. However, within a pMSSM scenario, there still 
exists some parameter space for light neutralino DM which could be completely probed by 
a dedicated analysis of the existing experimental data, in combination with the ongoing 
searches. 
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